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fNj ' Abstract 

> : 

\]q i We develop a unitarity approach to consider the final state interaction corrections 

Q\ • to the tree level graphs calculated from Chiral Perturbation Theory (\PT) allowing the 

inclusion of explicit resonance fields. The method is discussed considering the coupled 
channel pion and kaon vector form factors. These form factors are then matched with 
the one loop xPT results. A very good description of experimental data is accomplished 
for the vector form factors up to y/s < 1.2 GeV beyond which multiparticle states play 
a non negligible role. For the P-wave tttt phase shifts the agreement with data stands 
even higher up to y/s < 1.5 GeV. We also consider the isospin breaking effects due to 
the lo — p mixing as a perturbation to the previuos results. In addition the low and 
resonance energy regions are discussed in detail and for the former a comparison with 
one and two loop xPT is made showing a remarkable coincidence with the two loop 
xPT results. 
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1 Introduction 

The study of the pion vector form factor is an interesting problem mainly because pions are 
the lightest hadrons and hence they are common products in many experiments so that a 
good description of the pion electromagnetic form factor is often required. 

Many of the studies of this problem deal with some kind of modified vector meson 
dominance (VMD) Q when taking into account the effects of unitarity and final state 
interactions. This was done long time ago in ref. 0] and it was found that these effects 
show up not only as a modification of the bare p propagator but also as a change in its 
bare couplings. Other more phenomenological parameterizations are the ones in refs. ||, |4j, 
which basically account for the dressing of the bare p propagator and allow to add more 
resonances and parameters in order to fit the data up to high energies, hiding in this way 
extra effects as the presence of multiparticle channels, e.g. Att, ujtt etc... This latter kind of 
expressions are the ones commonly used in many experimental papers in order to fit their 
data and determine the resonance content. 
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It is well known that the low energy effective theory of QCD is Chiral Perturbation 
Theory {\PT) ||, |(| [jj]. Although this is a systematic way to express the QCD Green 
functions in terms of a power momentum expansion, unfortunately, it is valid only for 
low energies. Hence, if one is interested in higher energies nonperturbative schemes are 
unavoidable. Nevertheless, one should demand that, when used at low energies, these 
nonperturbative methods reproduce the low energy constraints of QCD given by the xPP 
expansion. 

The approach described in sections 2 and 3 reproduces the one loop xPP pion and kaon 
vector form factors Q and, as we show below, it is also appropriate to study higher energies 
by satisfying unitarity in coupled channels and incorporating explicit resonance fields |J. In 
this latter reference it is discussed how, at lowest order in the chiral expansion, resonances 
with spin < 1 couple with pseudoscalars (tt, K, rf) and with electroweak currents. It is 
also shown there how to connect with the standard VMD picture, when taking into account 
results from large N c and QCD high energy constraints. 

Other nonperturbative approaches that match with the one loop xPT vector pion form 
factor are given in refs. || 11]. In ref. [12, 13] the matching is given up to the two loop 
SU{2) x SU{2) xPP pion vector form factor first calculated numerically in [13] and then 
analytically in 14]. While the work of ref. [13] is only interested in assessing the relevance 
of higher order loops at low energies, the works of refs. 
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12 ] make predictions at 
Ref. 



higher energies, although all of them take into account only elastic unitarity. Ref. |9[, as 
1 13], makes use of an Omnes representation^ |la] but allowing for the explicit presence 
In ref. Ill]] the pion vector form factor is studied by solving the 



ref. 

of the p resonance 

Bethe-Salpeter equation making some assumptions on the off shell extrapolation of the 
amplitudes and using only the tttt channel. On the other hand, in ref. |l(], [l^] [0,1] and [0,2] 
Pade approximants from the xPT results are used (see also the discussions in refs. [^, [16| 
with respect the former references). Finally in refs. [17, [l8j an Omnes parameterization, 
assuming elastic unitarity, is implemented using directly the experimental 1=1 P-wave tttt 
phase shifts. 

In our work the isospin limit is taken but we shall also estimate the effect of the isospin 
violating uj — p mixing in the pion vector form factor. This contribution manifests itself in a 
very narrow energy range to the right of the peak of the p mass distribution in accordance 
with the mass and small width of the u (see fig.@). Other works in which this effect has been 
already discussed are [Q, [0^, pcfl . On the other hand, when comparing our results for the 
vector pion form factor in the p region with the experimental data, we shall also take data 
from T-decays [21] in which only the 1 = 1 part of the vector pion form factor contributes. 
Finally, in ref. p2] the isospin breaking effects in the pion form factors for the low energy 



regime are studied at next-to-leading order in SU(2) x SU(2) xPT- 



2 Tree level 

We will evaluate the tree level vector form factors and scattering amplitudes making use of 
the lowest order xPP Lagrangian plus the chiral resonance Lagrangians ||. These tree 
level amplitudes will be leading in large N c , while non leading contributions in this counting 
will be generated through the unitarization process to be discussed in the next section. 
In the case of the vector tttt and KK form factors we have to evaluate the diagrams 

1 Although in the end, since only the real part of the exponent of the exponential is kept, some of the 
analytical properties of the form factor are lost. 
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a) b) c) d) 

Figure 1: Diagrams considered for calculating the tree level form factors, eqs. ([|). From 
left to right: a) xPT lowest order, b) Octet, S, and singlet, Si, scalar resonance exchanges 
coupled to the vacuum, c) Wave function renormalization from the exchange of scalar 
resonances and d) exchange of vector resonances, p, to and (ft. Fig. |]a is calculated from 
ref. @ and the rest from ref. ||]. 



depicted in fig.Q. The contributions from the diagram of flg.|l]b and the ones from the wave 
function renormalization, fig|l]c, cancel each other due to charge conservation. A similar 
set of diagrams has been considered in ref. [24| in order to study the coupled Kir, Kr] and 
Krf scalar form factors [25]. In the former reference, although restricted to the study of 
the associated S-wave 1=1/2 meson-meson T-matrix, the requirement that the form factors 
vanish at s — > oo is used to reduce the number of free parameters. 

On the other hand, we will assume ideal mixing between the oj$ and u)\ resonances, so 
that the to and (ft are given by 



T^ 1 + 7T 8 



i 



--U)l 



v/3 1 



(1) 



and we will work explicitly in terms of the uj and (ft mesons. In this way we will be able 
to distinguish between the different physical masses of the u and (ft resonances. Note that 
at the order considered in ref. [|| vector singlet resonances, oj\, do not couple neither to 
photons nor to pseudoscalars. 

The resulting tree level vector form factors for the it + it~, K + K~ and K°K systems 

are: 
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(2) 



where Fy measures the strength of the photon-vector resonance vertex, Gy the same but 
for the pseudoscalar-vector resonance ones and s is the usual Mandelstam variable. Finally, 
Mp, M^ and M^ refer to the bare masses of the p, uj and (ft resonances and / is the pion 
decay constant, f w , in the chiral limit 0. 

One can proceed analogously for the evaluation of the tree level scattering amplitudes 
between the former states. In ref. [^] similar tree level amplitudes were already calculated 
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a) b) c) d) e) 



Figure 2: Diagrams considered for calculating the tree level amplitudes, eqs. (||). From 
left to right: a) xPT lowest order, b) Wave function renormalization from the exchange of 
scalar resonances, c) Tadpole-like diagram with the exchange of scalar resonances coupled 
to the vacuum, d) Exchange of vector resonances, p, lo and (j), in the s-channel and e) 
crossed exchange of vector and scalar resonances. Fig. ||a is calculated from ref. |?J and 
the rest from ref. Q. 



in the same way for studying the elastic tttt and Kir scattering. In addition, in ref. [24] 
the Kr] and Krf channels were also included to study the Kir scattering up to 2 GeV. 
The corresponding generic set of diagrams is indicated in fig.||. However, out of these 
diagrams we are not going to consider explicitly here those corresponding to the exchange 
of resonances in the crossed channels, fig]2]e. The reason is twofold: 1) they are not necessary 
to match with the one loop xPT results for the SU(3) vector form factors [f?j since they 
give rise to higher orders. 2) Because of VMD one expects a good description of the 
scattering amplitudes in the physical region when including the s-channel exchange of the 
vector resonances together with unitarity. With respect to the second point one has to 
take into account that in the chiral Lagrangians J?], || the standard picture of VMD |23|] 
is only accomplished when considering the lowest order xPT contribution together with 
the explicit resonance fields [^7], |9|. Furthermore, when making a dispersive analysis of the 
vector form factors, only the scattering amplitudes in the physical region are involved (see 
next section). 

The resulting P-wave partial wave amplitudes are: 
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By writing the lowest order amplitudes in terms of /, there is a cancelation between the 



"A partial wave amplitude with angular momentum L is defined to be 



dcos 6P L (cos 0)T(s, t, u) 



where Pl (cos 9) is the Legendre polynomial of Lth degree, T(s,t,u) is the scattering amplitude and s, t and 
u are the usual Mandelstam variables. 
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contributions of the wave function renormalization terms, fig. |2]b, those of the exchange 
of the scalar resonances coupled to the vacuum, fig. ||c, and the 0(p 4 ) crossed channel 
scalar contributions absorbed in the 0(p 4 ) xPT term proportional to the L5 counterterm 
J?], U, fig. ^e. Consistently with our approach of working the large N c leading contributions 
from ref. |7|], we have taken the masses of the scalar singlet and octet equal and we have 
also made use of the relation between their couplings deduced in ref. In this way, for 
instance, one has that L4 = Lq = consistently with phenomenology [28] and with the fact 
that both are subleading counterterms in large N c . 



3 Unitarization 

We will deduce our final expressions for the form factors after deriving those of the T- 
matrix. The latter will be accomplished by unitarizing the tree level scattering amplitudes 
in eqs. (|3|) following ref. |p7| . In this reference the most general structure of a partial wave 
amplitude when the unphysical cuts are neglected was deduced. It was also shown that this 
structure is the one required in order to match with the tree level amplitudes calculated 
from the lowest order xPT Lagrangian j/J plus the exchange of the resonances as in 
eqs. (g). In refs. [^], [^] the method is extended to include unphysical cuts as calculated 
in one loop xPT with and without baryons, respectively. For a review on these techniques, 
see ref. [|3l] . 

We work in the isospin limit. Our isospin states are: 



1 = 1 

1 !, + -- + 

\TTTT > = — \7T IT — 7T 7T > 

\KK > = —=\K+K- - K°K° > 

1 = 

\KK > = -^=\K + K- + K°K° > 



(4) 



note the extra factor in the normalization of the \irir > I = 1 vector since in this 

state pions behave like identical particles. We have also removed a global e l7T phase in all 
the states in order to have the form factors positive defined in s = 0. 

Since for I = 1 we have two coupled channels, we will use in the following a matrix 
notation in which the pions are labeled by 1 and kaons by 2 in the 1 = 1 case. For 1 = 
we have only the kaon channel which we denote by 1. 

Taking into account eqs. (||) and (f|) we can calculate for both isospin channels the 
tree level amplitudes between definite tttt and KK isospin states. We will collect these 
amplitudes, for each isospin channel, in a i^j-matrix. Thus, from ref. [27| we have the 
following expression for the T/-matrix: 

Tj(s) =[I + Kj(s) ■ g^s)] - 1 ■ Kj(s) (5) 

where I is the unity matrix and g^s) is a diagonal matrix [^7j given by the loop with two 
meson propagators. In dimensional regularization one has explicitly: 
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(6) 



where the sub index i refers to the corresponding two meson state and (Xi(s) = yl — 4m? /s 

with m, the mass of the particles in the state i. For the choice d\ = one has 5/ (s) = —J{s)u 
with the J(s)u functions defined in ref. such that J(0)jj = 0. 

We now introduce the diagonal matrix Q(s)a = Pi(s) 9(s—4mf), withpj(s) = \J ' s/A — mf 
the modulus of the cm. three momentum of the state i and 9{x) the usual Heaviside func- 
tion. In terms of this matrix we have 

Sj(s) = I + » Q(s) 1 ' 2 ■ Tj(s) • Q{s) 1 ' 2 (7) 

47TV S 

with Sj(s) the S-matrix fulfilling the requirement Sj(s) ■ Sj(s)^ = I. 

The electromagnetic meson form factor, i*MM'( s )> ^ s introduced as follows. The transi- 
tion amplitude from a photon, virtual or real, to a pair of mesons is written as: 

< 1 {q)\T\M(p)M'{p') >= ee^p-p'TFMM^s) (8) 

with q 2 = s, e the modulus of the charge of an electron and the photon polarization 
vector. 

Unitarity of the 5-matrix expressed between the 7 and the MM' states, up to C(e 2 ), 
leads to the relationship: 

Pmw(s) ImF(s) MM > = F a (s)^jLp a (sMs - Am 2 a ) (T^V(s))* (9) 

a » 

where * means complex conjugation and the strong amplitudes are projected in the P- 
wave. In the following we will remove the supraindex L = 1 with the understanding that 
any strong amplitude will be in P-wave. 

Dividing the former equation by Pmm'( s ) an d taking the complex conjugation of the 
term in the right hand side, since this is a real quantity above the threshold of the system 
\MM' > where unitarity applies, one has: 



ImF MM ,(s) = J2F* a (s)^^e( S - Am 2 a )p a (s 



a, MM' 



PMM'{s) 



(10) 



Taking into account that the T(s) a p is a symmetric matrix with respect indexes a and f3 
because of time reversal invariance, introducing also the diagonal matrix Q(s)ij = Pi(s) 5ij 
and the one column matrix F(s)i = Fi(s), we can write eq. ( |Iol) in a matrix notation as: 



ImFOO = Q(s)- 1 ■ T(s) ■ • Q(s) ■ F*(s) (11) 

If we further substitute in the previous equation ImF(s) by (F(s) — F*(s))/(2i) and 
T(s) by its expression in eq. (||), one finds: 

F(s) = ^I + Q( s )-^[I + K(s)-g(s)]- 1 -K(s)-z^-Q(s)YF*( S ) 
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= Q(s)- 1 [I + K(s) • g(s)]- 1 • { [I + K(s) • g(s)] ■ Q(s) + K{s) ■ t^j= ■ Q(s) } • F* (s) 

= Q(s)- 1 ■ [I + K(s) ■ g(s)r l ■ Q(s) • {/ + Q(s)" 1 ■ K(s) ■ Q(s) ■ g(s) 

+ Q( s r 1 -K(s)-Q(s)-i^=}-F*(s) (12) 

In the last equality we have made use of the fact that the matrices Q(s), Q(s) and g(s) 
commute since all of them are diagonal. 

Let us a first note that the symmetric matrix K(s), eq. (|3|), is a matrix of functions with 
a kinematical cut in its K(s)i2 matrix element between the threshold of channels 1 and 2 
because K{s)\2 is proportional to P\(s)p2{s). However, the matrix Q(s)~ 1 • K(s) ■ Q(s) has 
no such kinematical cuts since the modulus of any of the three-momenta p a appears always 
squared. Note also that, although does not exist for s equal to its value for any of 

the thresholds, by continuity, the product Q(s)" 1 • K(s) ■ Q(s) is well defined, and hence 
Q(s)^ 1 • K(s) ■ Q(s) is a well defined real matrix. Taking this into account and that 

9{sr=9{s)+l ^=s (13) 

since 

W») = -||$, (14) 

we can write eq. (|l^) as: 



F(s)= I + Q(s)- l -K(s)-Q(s)-g( S ) ■ I + Q(s)- 1 ■ K(s) ■ Q(s) ■ g* (s) -F*(s) (15) 



Multiplying both sides of the former equation by 
arrive at the interesting result: 



wc 



I + Q( s )-l.K( S )-Q(s)-g(s 
I + Qisy 1 ■ K(s) ■ Q(s) ■ g(s) ■ F(s) = I + Q(s)- 1 ■ K(s) ■ Q(s) ■ g*(s)] ■ F*(s) (16) 



This implies that the matrix 



F(s) has no cuts at all, since 



I + Q(s)- 1 -K(s)-Q(s)-g(s) 

the right hand cut or unitarity cut, the only one present in g(s) or F(s), has been removed. 
Note that, as it is well known, a form factor is an analytic function of s except for the 
presence of the right hand cut from threshold up to infinity. It can also have poles on 
the real s axis below threshold, corresponding to bound states, and poles in the complex 
plane in the unphysical Riemann sheets corresponding to resonances. On the other hand, 
as discussed above in contrast to K(s), Q(s)^ 1 ■ K(s) ■ Q(s) is a function free of physical or 
kinematical cuts. Prom eq. (|l~6| ) we write now 



F(s) 



I + Q(s)- 1 -K(s)-Q(s)-g(s)\ ■ R(s) 
Q(s)- 1 -[I + K(s)-g(s)r 1 -Q(s)-R(s) 



(17) 



with R(s) a function free of any cut. Note also that the matrix D(s) = I + K(s) ■ g(s) was 
already introduced in eq. (||) in relation with the purely strong T-matrix. 
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As told above, the considered large N c leading contributions are given in eqs. (§). In 
fact in this limit D(s) = I + 0(N ( T 1 ) and hence we have: 

F( S ) = R( S ) Nc l eading = F* (s) (18) 

where F t (s) is the tree level form factor. As a result we can write: 

R(s) = F t (s) + SR(s) (19) 

with 5R(s) subleading in large N c and free of any cut. Since the pole singularities present in 
the form factor correspond to tree level resonances (already present in F t (s)) or dynamically 
generated resonances from the D(s) matrix, see eqs. ([l?]) and (|5|), the SR(s) matrix should 
then be simply a matrix made up of polynomials^]. 
Hence, we can rewrite eq. (|l~7|) as: 

F(s) = \l + Qis)" 1 ■ K(s) ■ Q(s) ■ g(s)} • [F\s) + 5R(s)] (20) 



We can further reduce eq. ( |20|) for the two channel case. Indeed, given the structure of 
the K matrix of eq (0) for 1=1 and -kit and KK channels we find 



FM _ 1 ( 1 + P2(s) 2 0(s) gl(s) -V2p 2 (s) 2 f3(s)g 1 2 (s)\ ( t , 
HS) ~ det[D(s)] [ -V2 Pl {sf 13(8) g\(s) 1 + 2 Pl (sf 0(b) g\( 8 ) J ^ W + dH W) ^ 

with (5 given by 

m = £ (i + 2 fj^) m 

and det[D(s)] = 1 + 2pi(s) 2 f3(s) g{(s) +P2(s) 2 f3(s) g^is)- For 1 = 0, with only one channel, 
eq. (|2(i| ) is simple enough. In the explicit formula of eq. (|2l| ) we can see that there are 
indeed no kinematical cuts associated to the P1P2 factors of the matrix K, as argued above. 

The unitarity method we have introduced is similar to the one discussed in ref. |32|] 
in order to solve the Muskhelishvili-Omnes problem [15]. The differences arise from the 



fact that we have showed how to deal with the kinematical singularities associated with 
the product pi"hnF(s)i eq. (||), and also because we have allowed the presence of pole 
singularities both in D(s), due to K(s), and in F(s). However, they appear in such a way 
that they cancel each other and the results are finite. In that reference [^2| the connection 
between eq. ( |20| ) and an Omnes representation for the elastic case is also discussed. This 
relation is easily seen since in the elastic case D~ 1 (s) has the phase of the strong scattering 
amplitude, eq. (|5]), and only has the right hand cut. There is, however, a subtle point which 
is that D~ 1 (s) has zeros for s equal to the bare masses of the vector resonances introduced 
in K(s) and they have to be removed in order to make an Omnes representation. This is 



accomplished automatically in eq. (20) due to the explicit presence of the bare resonance 
poles in F l (s). 

The degree of the polynomials present in 5R(s) can be restricted when considering the 
behavior of the form factors in the limit s — ► oo. In fact from eqs. (|2|), (||) and ( |20| ) we have 
the following limits: 

3 It should be understood in the former sentence that subleading corrections in 1/N C to the couplings and 
masses of the resonances (not coming in our formalism by the D(s) function), have to be reabsorbed in a 
renormalization of F t (s), eq. (H). 
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Figure 3: W is denned as ^fs for s > and as —y/—s for s < 0. From left to right: a) ir + ir~ 
vector form factor. Experimental data from ref. Q collecting also data from refs. p8[ for 
positive s. For negatives s ref. [39] b) -kit P-wave phase shifts. Triangles from ref. |40| and 



circles from ref. [41| 
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2G? 
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(23) 



where g\ (00) is a quantity logarithmically divergent with s — > 00, see eq. (|6|). We denote by 
(00) the value of the vector element of 5R(s) with isospin / and channel i for s —* 00. 
Clearly, this value diverges as s n<yI '^ with n(I;i) the degree of the polynomial 5Rj (00). 



From perturbative QCD 33], it is known that the vector pion form factor goes at most 
as a constant for s — > 00. Although experimentally, and also from the quark counting 
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rules, it is very likely that it vanishes for s — > oo. Experimentally, we have Fy = 154 MeV 
from the observed decay rate r(p° — > e + e") || and Gy = 53 MeV from a study of the 
pion electromagnetic radii after taking into account corrections from chiral loops at next 
to leading order With these values, we see from eqs. (23) for the 1 = case, that 
5Ri(s) can have, at most, degree one since otherwise the form factor would be divergent 
(we are assuming here that the high energy limit of the kaon form factors is the same as 
the aforementioned one of the pion vector form factor). Analogously for / = 1, we see that 
the difference 6R{(s) - V25R\(s) can be at most a constant, while, independently, 5R\(s) 
can have degree one. If we further require that the vector form factors vanish for s —* oo, 
both for the pion and kaon ones, then 5R\ (s) are constants and 



5R\ = V25R\. 



(24) 



In the following we will assume that this is the caseQ. Note that since 5Rj(s) are at least 
0(p 2 ) quantities, these constants must be proportional to the square of the pseudoscalar 
masses m 2 , m 2 K . Then we will finally have: 



F(s) = \l + Q(s)- 1 • K(s) • Q(s) • g(s)} • (F*(s) + a T ) 



(25) 



where ai is just a constant for 1 = and for 1 = 1, taking into account eq. (24), a\ can be 
written as: 



O'l 



a 



'iA/2 



(26) 



with cti a constant. 



3.1 Matching with one loop x_PT 

In order to fix the values of the constants d\ of g\[s) and aj of eq. (25) we match our 
results with those of xPT up to one loop ||. The form factors from together with eq. 
(H) for the definition of the isospin states, are: 



-^ivr 1 = 1 + 2Hwn + HKK 



7T7V 
7,1=1 



Fj^ = -±=(1 + 2H^ + H KK ) (27) 



1 

3 V 



F K = R = -^(I + ZHkk) 



with the function Ha(s) given by: 



Hu(s) - t^L'M + ^ (2^,? - atl.) log ^|±i - 1 + log £,) (28) 

In order to make the matching with the one loop xPT form factors, we will take into 
account the fact that for evaluating the contribution of the octet of vector resonances to 
the 0(p 4 ) xPT counterterms, for instance to Lg(/x), the masses of these resonances have to 

4 We refer the reader to refs. ^) which stress the importance to keep the short distance QCD constraints 
in order to guarantee the independence of the results under spin-1 field redefinitions. 



10 



be taken in the chiral limit ||. We will denote this mass in the following by My ~ M p ||. 
As noted in ref. Q the Lg(p) is saturated by the meson vector resonance exchanges at a 
scale fi around the mass of the p, such that 

2L 9 (p « Mfy* ical ) = (29) 

we will also make use of this result in the following. 

Let us consider first the 1 = 1 case. Expanding our results, eq. (25), up to one loop in 
the xPT counting, we have: 

?I=lr„\ i , , F vG v s s f 2 2 , 3 <7 w (s) + 1 



F<->) = l + ai + ^ + _| 2 a,(^- < r,(.)^-M.) J lcg^- I | + 

+ w{ 2 -W 2 --W^--W 3 '^S|5fI} < 30 > 

while the one loop X-P^ result is: 

^ lSj " i+ /2 + 96tt 2 / 2 1 g ^( S )-l 3 ° g /i 2 J + 

+ T ^{2^) 2 -^) 3 lo g ^4^-i-log4| (3D 
192tt z j z { cfk(s) — 1 3 fi z J 

The matching for KK form factor in / = 1 does not give any new condition since this 
form factor, at the one loop chiral level, is simply l/v2 of the 7T7t one, both in our approach 
and in xPT ]?]]. Note that this happens independently of the value of the constant ctq, 
because of eq. fl26|). 

In the following we will take 

a[ = (32) 

in order to constraint further our approach. This can be done since, as we show below, we 
can match our results with one loop xPT by choosing appropriate values for d\ and d\. On 
the other hand, as we will see in the next section, a very nice description of the pion form 
factor and of the ttit P-wave phase shifts, see fig. ||, is given without including this extra 
degree of freedom. In addition, we will discuss below the sensitivity of our results under 
changes of ol x ^ 0. 

By matching eqs. ( j30l ) and ([31]) with a\ = 0, taking also into account eq. (p9|), one 
finds the condition: 



^{s?d\ + °-^ftd\ = \ + log 4 + \ 4 (33) 
Identifying the terms independent of s and linear in 1/s, we finally have: 



^^^flog^ + ^og^- 1 



m? K — m 2 \ /i 2 2 /i 2 2 



,i -2 m 2 / ml 1, m\ 1\ 



mi — 7774 \ /i z 2 p z 2 
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Next we do the matching in the 1 = sector. Our chiral one loop expression for this 
form factor is: 



V2Fj=? = 1 + V2a + + S 



M 2 / 2 64vr 2 / 2 



The xPT one is: 



2a K (sf - d\a K {sf - o K {sf log + 1 1 

cr^(s) -1J 

(35) 



By matching eqs. 
condition: 



2a K {sf - a K {sf log + \ - ~ - log ^ }> (36) 



1 



and (|36|), taking also into account eq. (29), we obtain the 



\/2a 



sax 



{sfd\ 



1 



/?? 



K 



64vr 2 / 2 V3 +1 ° g V 2 



64^2/2 

Identifying the terms independent of s and those linear in s: 



(37) 



m 



a 



K 



o + lo g 



77? 



A' 



1 3 ' ^ 6 ^ 



16^2 vr 2 / 2 V3 



1 . mi 
+ log^ 
A* 



(38) 



We see indeed that ao is subleading in N c since f 2 ~ A^ c and m? K ~ 1 in this counting. 
Numerically we find it to be of the order of 0.1 for \i = M p , consistent with the large iV c 
counting. 



4 Results 

Once we have fixed the constants d\ and aj by eqs. (]32]), (|J) and (|3g|) , our final expression, 
eq. (25), for the vector form factors only depends on the bare masses of the p, cj) and u 
resonances and on the couplings Gy and Fy and on /. The couplings Fy and Gy are fixed, 
as explained in section 3, from their experimental values, 154 MeV and 53 MeV, respectively. 
The value for the parameter / is taken from the second entry of ref. j/J , where it is derived 
by working out the relation of the ratio f n / f with the isoscalar uu + dd radius of the pion. 
The resulting central value is / = 87.4 MeV. Notice that this estimate is done for the value 
of / in the limit m u = rrid = and m s ^ and there is some controversy about the possible 

. The bare masses 



deviations with respect this value when considering the SU(3) case [35 



can be fixed in terms of Fy, Gy and / by the requirements that the moduli of the ttit 1 = 1 



and KK 1=0 P-wave amplitudes have a maximum for y/s = Mj 



physical 



770 MeV and 



for ^/i = Mf ysical = 1019.413 MeV, respectively. We obtain the values M p = 829.8 MeV 
and M,/, = 1026.581 MeV, where the number of decimals correspond to the experimental 
precision in which the physical masses are given [37]. 

For the mass of the u we take directly 782 MeV since there are no experimental data in 
the region of the u and its contributions to other physical regions do not depend on such 
fine details since the uj is very narrow. 
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Finally, note that our physical results does not depend on the regularization scale \x 
since any change in the scale is reabsorbed by the corresponding change in the constants 
df(/f)- However, since we have made use of eq. (|29|), we have to consider a value for p, 
around the mass of the p meson, where vector meson saturation of the L^ip) coefficients 
works |Q. In the following we will fix p = M p . 



3 - 




500 1000 1500 2000 1000 

W(MeV) 



1200 



1400 
W(MeV) 



Figure 4: W is defined as ^/s for s > and as — \/— s for s < 0. From left to right: a) 
K + K~ vector form factor. Experimental data: black circles fl42j , white circles Q and stars 
|| . b) a(e + e~ K°K°) in nb. Data: g| 

As can be seen in figs. [3], |] and [|, we can describe in a very precise way the vector 
pion and kaon form factors from negative values of s up to about 1.44 GeV 2 . The P- 
wave 7T7T phase shifts are also very well reproduced even for energies higher than those for 
the form factors, fig. || is displayed until ^fs = 1.5 GeV. Note that the curves present 
in those figures do not result from a fit making use of our approach since the parameters 
have been fixed in advance. Furthermore, the qualitative behavior of the data for the form 
factors is reproduced even for very high energies. For values of -^s higher than 1.2 GeV 
new effects appear: 1) the presence of more massive resonances, p' , oj' , <p' ... These states 
can be in principle taken into account by our formalism just by adding more resonance to 
eqs. (|2|) and (||). 2) More seriously and less trivial is taking into account the effect due to 
multiparticle states, e.g. 47r, u>tt ... In fact, since in our model the widths of the resonances 
are derived dynamically in terms of the included channels we cannot mimic the effect of 
such multiparticle states as done in other studies in which the widths of the resonances are 
fitted, that is, incorporated by hand. This is currently done when using parameterizations 
of the type of ref. [||, ||]. In fact, in the present case we have found that the effect of other 
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channels cannot be neglected for energies higher than ^/s > 1.2 GeV, for which we do not 
give a fair reproduction of all the structures present in the data. We have checked that the 
inclusion of higher mass resonances like the p' improves the agreement with the data but 
there are clear signs that other elements are still missing. 

In fig^a one clearly sees the peaks of the 4>, p and u resonances, the latter on the top 
of the p. The uj peak corresponds to that of a zero width resonance since it appears below 
the KK threshold and the 3tt state is not included. 

In figs.[|a,b we compare for low energies our results for the pion and charged kaons vector 
form factors with the results obtained in x_PT for the pion vector form factor, up to one |?J 
and two loops fl3| , |l4|], and with the one loop xPP |?j] charged kaon vector form factor. In 
the figures the aforementioned matching, guaranteed by construction, is clearly seen, and 
the resummation accomplished by eq. (p5|) provides the right corrections to the xPT results 
also for low energies. In fact, we see a much better agreement of our result with the two 
loop xPT pion vector form factor than with the one loop result. 

In table 1 we give the calculated electromagnetic radii of the pions, charged and neutral 
kaons and compare them with one loop xPP J/J and with experiment. We have not shown 
the two loop xPT value for < r 2 > since its experimental value [48] is taken as an input in 
order to fix a counterterm [13] . In ref. [14] it is argued that this is a sensible assumption 
within an accuracy around a 10% in the value of < r 2 >y. Our results for pions and 
charged kaons are compatible with the xPP results and experiment within errors. For the 
case of neutral kaons, assuming as we have, the same uncertainties as in the xPT results, 
we are also compatible with experiment. It is also interesting to compare our values for 
the electromagnetic radii with a low energy theorem due to Sirlin [ 51 1 . This theorem states 
that, 

1 

2 

where m s is the mass of the strange quark, fa is the average between the masses of the up and 
down quarks and f+ n is one piece of the vector Ktt form factor j/J . In the previous relation 
the contribution of heavy quarks have been neglected and we follow the notation of ref. @ . 
Numerically at 0(p 4 ) in xPP Rj 3116 nas ' f° r ^ ne same combination of electromagnetic radii 
as in the right hand side of eq. @), 0.00±0.05 fm 2 . At 0{p 6 ) the result is |2| 0.021 ±0.003 
fm 2 . From the values shown in table 1, taking from experiment the electromagnetic radius 
< r 2 >f 7r = 0.36 ± 0.02 fm 2 |7), as also done in 0(p A ) xPT, we have -0.07 ± 0.05 fm 2 . 
The experimental value, from refs. f|| |5(], |7||, is -0.02 ± 0.05 fm 2 . 



0[(m s - fa) 2 } = ± (Fy(s) + F« + (,)) + F*\s) - (39) 



In figs.||c,d the resonance regions for the p and (f> mesons are shown in detail. In the case 
of the pion form factor in the p region, we have compared our results with the experimental 
data of ref. [21] since in r decays only the 1 = 1 part of the pion form factors plays a role 
and we are working in the isospin limit. As can be seen the reproduction of data is very 
good in both cases. 

As announced before, we have also studied the dependence of the results on the sub- 
leading constant a' x which has been set to zero in eq. p2]). Notice that the matching to 
xPT gives new values for the d\ parameters in terms of those of and also one has to 
recalculate the bare mass of the p resonance. The new results, obtained with reasonable 
values of q' 1; can be hardly distinguished from the a\ = case in the region y/s < 1.2 
GeV. A fit to the tttt phase shifts and the form factors data gives that the best values are 
a'i = —0.0012 and M p = 830.23 MeV, leading to results that are indistinguishable from the 
ones obtained by setting a[ = 0. 

Finally, in our approach we have taken into account tttt and KK coupled channels for 
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the 1=1 P-wave. The effect of coupled channels is more important for kaons, particularly 
for the neutral ones, than for pions. This can be expected since when we decouple the 
kaons from the pions the p resonance appears in the kaon vector form factors as a zero 
width resonance and this is a very bad approximation, particularly, if we are interested 
around the rather broad p energy region. On the other hand, around the peak of the (ft 
the effect of the pions is small, of the order of a few percent, and the kaon vector form 
factor is dominated by the KK elastic channel strongly coupled to the (ft. However, as we 
move again to higher energies, where the form factors have dropped substantially, the pion 
channels become again important. 




500 600 700 800 900 1000 1010 1020 1030 1040 

W(MeV) W(MeV) 



Figure 5: TV is defined as yfs for s > and as — \J— s for s < 0. From left to right and top 
to bottom: a) Vector pion form factor. The dashed-dotted line represents one loop xPT 
ref. [J/J and the dashed one the two loop xPT result ref. b) K + K~ form factor. The 
meaning of the lines is the same as before. Data are from ref. ^6|, white circles, and the 
star points are from ref. ||47||. c) Vector pion form factor in the p region without taking into 



account the p — u> mixing effect. Data from r-decays |21[|. d) K + K form factor. Data [42] 



4.1 Isospin violation: p — u mixing. 

The incorporation of isospin violation effects, due to the mass differences from the u and d 
quarks and electromagnetic corrections, is readily possible, up to some extend, within the 
present formalism. As an example, consider the OZI violating process (ft vr + 7r~ studied 
in ref. [53| within a chiral unitary approach. It was noted there that the u — > tt + it~ 
decay requires a direct p — u> mixing which makes this process qualitatively different to 
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e.m. radii (fm 2 ) 


Our results 


One loop X-PT 


Experiment 


< r 2 >y 


0.44 


0.44 ± 0.04 


0.439 ± 0.030 [50] 


< r 2 


0.34 


0.38 ± 0.03 


0.34 ± 0.05 fT| 
0.28 ± 0.07 @ 


< r 2 


-0.10 


-0.04 ± 0.03 


-0.054 ± 0.026 H 



Table 1: Electromagnetic radii for 7r + 7r - , K + K~ and i^°iir , from left to right respectively. 
In the table are given our results, the ones from xPT at next-to- leading order ref. |7| and 
some experimental data. The uncertainty in our values has to be considered of the same 
size than the one in xPT, that is, around ±0.04 fm 2 . 



7T + 7r . This p — uo mixing has been studied in the context of chiral Lagrangians, 



the c, 

together with the large N c expansion, in 54] with the result: 



with 



e 



pui 



(m 2 K0 



m 2 K +) + (mlo 



(40) 



+ 



e 2 F 2 



(41) 



Given the smallness of the contribution from p — lo mixing one can add to the result for 
the pion form factor evaluated before the contribution from the mechanism of fig.|6[ 

co p / 

r \T " x c 



Figure 6: Feynman diagram corresponding to the p — co mixing contribution to the pion 
form factor. 



The contribution of fig|] is straightforward and has been evaluated previously in [20] 
With the ■joj and pinr couplings from || and the p — lo coupling of eq. (|40|) we obtain: 



AF 



I FyGy 1 



s-Mj + iy/sT p (s) 







(42) 



Since this correction is only visible around the p peak we take for T w (s), Tp(s) the constant 
values of the PDG [[?<]] (using instead an energy dependence width of the p meson does not 
introduce any appreciable effects). 

With a value for Q pui of —4500 MeV 2 from eq. (|4l"|), which is also consistent with 
empirical determinations [55, 56], we get the form factor shown in fig]^ which reproduces 
fairly well the experimental data at the peak of the pion form factor. It is worth mentioning 
that this effect of isospin violation is not observed in the data from r decay since this process 
is only sensitive to the 1=1 current (see figff). 
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Figure 7: Solid line: pion vector form factor taking into account the p — uj mixing; dashed 
line: 1=1 vector form factor in the isospin limit. Data: black circles from e + e~ — ► 7r + 7r - 
reactions @, 38]; white circles: data from r-decays ]^j}] . 



5 Conclusions 



We have developed an appropriate method to take into account the final state interaction 
corrections to the tree level amplitudes of refs. M, |8| with explicit resonance fields. 

The strong tree level scattering amplitudes from refs. [Q, § were unitarized in coupled 
channels following the scheme developed in ref. ]p?|] . These amplitudes were then imple- 
mented in order to derive the final state corrections to the tree level vector form factors, 
eqs. (||). As a result a very good description of data is accomplished for yfs < 1.2 GeV. 
For higher energies multiparticle channels are no longer negligible and also new resonances 
appear, although our results still give the trend of the experimental data. 

Our expressions reproduce the xPT expansion of the pion and kaon vector form factors 
up to one loop Q. The resummation of our scheme leads, in the low energy region where 
xPT is expected to hold (see fig.|5|a), to a much better agreement with the two loop xPT 
pion vector form factor than with the one loop one. However, at the same time, we are 
able to provide a very good description at higher energies, including the resonance regions 
where standard xPT could not be applied, figs. ||c and |]d. Furthermore, our calculated 
electromagnetic radius of the 7r + , K + and K°, table 1, are in agreement with experiment, 
assuming the same uncertainties as in the X-PT results. 

We have also taken into account the p — u mixing as a main source of isospin breaking 
in order to compare our results with the form factor obtained from e + e~ — > 7r + ir~. The 
input for the p — uj mixing has been taken from recent work respecting also chiral symmetry 
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in the same line as the rest of the Lagrangians employed here. 

Finally, our method also provides P-wave scattering amplitudes for pions and kaons in 
very good agreement with the experimental data on the phase shifts for the former one. 
For the case of the kaons there are no scattering data available, but they would be highly 
desirable as a further test of the chiral unitary approach followed here. 
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